ON THE EXISTENCE AND UNIQUENESS OF WEAK SOLUTIONS 
FOR A VORTICITY SEEDING MODEL 
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Abstract. In this paper we study the Navier-Stokes equations with a Navier-type boundary 
condition that has been proposed as an alternative to common near wall models. The boundary 
condition we study, involving a linear relation between the tangential part of the velocity and the 
tangential part of the Cauchy stress- vector, is related to the vorticity seeding model introduced in the 
computational approach to turbulent flows. The presence of a point- wise non vanishing normal flux 
may be considered as a tool to avoid the use of phenomenological near wall models, in the boundary 
layer region. Furthermore, the analysis of the problem is suggested by recent advances in the study 
of Large Eddy Simulation. 

In the two dimensional case we prove existence and uniqueness of weak solutions, by using rather 
elementary tools, hopefully understandable also by applied people working on turbulent flows. The 
asymptotic behaviour of the solution, with respect to the averaging radius S, is also studied. In 
particular, we prove convergence of the solutions toward the corresponding solutions of the Navier- 
Stokes equations with the usual no-slip boundary conditions, as the small parameter 6 goes to zero. 

Key words. Navier-Stokes equations, boundary models for turbulent flows, existence, unique- 
ness, LES models. 
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1. Introduction. In this paper we consider the Navier-Stokes equations and in 
particular the role of boundary conditions in the simulation of boundary effects in 
turbulent flows. We consider the Navier-Stokes equations (in non-dimensional form) 
for viscous incompressible fluids in a bounded smooth domain Vl C R", n = 2,3, 

2 

d t u V-B(u) + (it- V)u + Vp = inftx(0,T) 

Re (i.l) 

V • u = in fi X (0,T). 

We recall that u — (u\, . . . ,u n ) is the unknown velocity field, p is the hydrostatic 
pressure, Re > is the Reynolds number, and E)(u) is the deformation tensor, i.e., 
the symmetric part of the matrix of derivatives of u: 

m , , If dui duk 

and the Navier-Stokes equations are generally equipped with the no-slip boundary 
conditions on T = dil 



u = Q onrx(0,T). (1.2) 

To introduce the problem that we will study, we recall that while at a free surface it 
is natural to require continuity of the stress-tensor (I being the identity) 

T(u,p) = -pI+-J-B(u), 
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the conditions at a solid wall are much troublesome. The no-slip condition 1|1.2|) has 
been justified by Stokes 30] since the contrary assumption 

. . . implies an infinitely greater resistance to the sliding of one portion 
of fluid past another than to the sliding of fluid over a solid. 
It is well-known that there are situation in which the boundary condition (|1.2|l may 
not be valid. For instance in Scrrin |28| §64 it is pointed out that in high altitude 
aerodynamics, or more generally when moderate pressure and low surface stresses 
are involved, the adherence condition is no longer true, see also the review in Trues- 
dell | 31) . In this respect several authors proposed various slip (generally nonlinear) 
conditions, modeling precise physical situations; see for instance Serrin [28], Beavers 
and Joseph Krein and Laptev [T8|. 

From the historical point of view, the slip (with friction) boundary condition 
proposed 1 by Navier was 

u-n = and /? u T + T(u,p) = 0, > 0, onrx(0,T), (1.3) 

where n denotes the unit normal vector to T, u T = u — (u ■ n)n, while T_( u tP) — 
t{u,p) — (t(u,p) -n) n denotes the tangential part of the Cauchy stress vector t defined 
by: 

n 

t(u,p) =n-T(u,p) = ^T lfc (u,p)n fc . 

k=i 

The controversial between the condition (|1.3|) proposed by Navier, versus the Dirich- 
let (|1.2fl proposed by Stokes was analysed also by Maxwell, who observed that the 
same conditions may be derived in the kinetic theory of gases and that the parameter 
(3 should depend on the Reynolds number Re and on the mean free-path A, satisfying 
the couple of consistency conditions 

[3 — > oo as A — > 0, for Re fixed 
[3 — ► as Re — > oo, for A fixed. 

With the above asymptotics it possible to recover in both cases the correct no-slip 
boundary condition for viscous fluids and the no-penetration condition for ideal fluids. 
A study of the numerical problems related with the implementation of (|1.3|l can be 
found in John 

Recently, Fujita ^I] performed the analysis with the "slip or leak with friction" 
boundary conditions. These conditions are of particular interest in the study of poly- 
mers, blood flow, and flow through filters. The boundary conditions studied in 
with the techniques of variational inequalities, turn out to be particular cases of the 
nonlinear boundary condition proposed at p. 240 in [25] and they are very strictly re- 
lated to both the Navier and the no-slip boundary conditions. See also Consiglieri jS] 
for similar problems. 

Among other nonstandard boundary conditions we recall that studied by Begue et 
al. @] and the "do-nothing" Neumann conditions, very appealing for numerical studies, 
implemented in Heywood et al. 

1 Note that in contrast to Stokes (1845) that used the continuum mechanics, Navier (1823) derived 
the equations by using some formal analogy with the elasticity theory and the assumption that 
molecules are animated by attractive and repulsive forces. 
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1.1. Near Wall Models and turbulent flows. Our interest in non standard 
boundary conditions comes essentially from the study of turbulent flows. First, we re- 
call that in the boundary layer theory several log-law and power-law are introduced, 
together with the fictitious boundaries, in order to model turbulent flows within a 
small region near the boundary. Roughly speaking, appropriate nonlinear boundary 
conditions are imposed on an artificial boundary that lies inside the computational do- 
main. The boundary conditions may simulate (at least in a computational approach) 
the behaviour of the boundary layer and they are modeled to take into account the 
peculiar behaviour of a fluid near the boundaries 2 

Our main motivation comes from the mathematical theory of Large Eddy Simu- 
lation (briefly, LES). In fact, the purpose of LES is to model the evolution of large 
coherent structures (eddies); this is done by studying the equations satisfied by a 
filtered velocity. Generally, the filtered velocity u is defined through a convolution 



with a rapidly decreasing smoothing-kernel g$ of width 5; in several cases of practical 
interest 



By definition, the value of u at a point xq on the boundary L will depend on the 
behaviour of u in a neighborhood of width 5 near that point; even if u is extended to 
zero for each x ^ fi, it is clear that in general u{xq) ^ 0. 

As pointed out in Galdi and Layton ^3] the physical intuition may suggest that 
large coherent structures touching a wall do not penetrate, sliding along it, and losing 
their energy. The boundary condition of Navier may be revisited, by linking the 
micro-scale A of the kinetic theory of gases, with the radius S of the averaging filter. 

Many Near Wall Models (NWM) have been tested in the computational approach, 
see Sagaut |2Z] and Piomelli and Balaras • The results are not uniformly successful 
and a positive application is very often based on a fine tuning of parameters. This is 
why new models require at least a positive background from the physical hypotheses 
and a coherent mathematical analysis. In particular, a successful application of the 
Navier slip-with-friction boundary condition i|1.3[) is prevented by two main facts: 
1) the presence of recirculation regions and 2) the presence of fast time-fluctuating 
quantities. 

The first problem is motivated by the fact that in recirculation regions the local 
Reynolds number is very different from the main stream and it is natural to expect 
that (3 should depend (possibly in a nonlinear way) on a local Reynolds number related 
to the local slip speed, i.e., 



Preliminary analysis has been performed by John et al. |17j . Dunca et al. |10 |. and 
an appropriate power-law choice of j3 seems promising to improve the estimation of 
reattachment points. 

The limitation of Navier law H1.3|l in a Boundary Layer theory is that it can well 
describe time-averaged flow profiles, while the information coming from fluctuating 

2 In this respect we quote J.C. Maxwell: "It is almost certain the the stratum of gas next to a 
solid body is in a very different state from the rest of the gas" . 



u{x, t) = gs(x) * u(x, t) 




0=m ki). 
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quantities in the wall-normal direction can play an important role in in triggering sep- 
aration and detachment. To try to overcome this limitation, very recently Layton [T§| 
recognized a particular class of boundary conditions, leading to conditions similar "in 
spirit" to the so called vorticity seeding methods. In fact, a Navier slip-with-friction 
boundary condition implies the generation of vorticity at the boundary, proportional 
to the tangential velocity. More precisely, in the case of a two-dimensional domain fi, 
for each smooth function v such that v ■ n — on the boundary, it holds 

n • D(v) ■ t — —curl?; + k(v ■ r) = on T, 

where curlu = dvi/dx 2 — dv^jdxx, r denotes the unit tangent vector, while k is the 
curvature of T. 

In particular, in |19| the following boundary condition is proposed to simulate the 
boundary effects 

u ■ n = 8 2 g(x, t) and - — u r +T(u,p) = 0, (1-4) 
oRe 

where g is a highly oscillating function in the time variable (hopefully a random 
variable in numerical tests), while it may be very smooth in the space variables and 
it should satisfy the natural compatibility condition 

g(x,t)da = Vie(0,T), (1.5) 

r 

which is required by the normal trace of a divergence-free vector field. 

This way of reasoning is also similar to the introduction of stochastic fluctuations, 
to simulate the micro-scale effects. A comprehensive introduction to stochastic partial 
differential equations in fluid mechanics and the statistical approach can be found in 
Monin and Yaglom |21| and one main mathematical paradigm is that an additional 
non-smooth term on the right-hand side may naturally take into account the effect of 
the fast fluctuating quantities. This leads to study the system 

d t u - J-V • D(«) + (u • V)u + Vj> = / + 

Re at 

where g is a function that does not have a proper time derivative, but is just continuous 
or with other weak properties. Appropriate notion of solutions, together with the 
statistical properties of the above problem are studied in Bensoussan and Temam [B] 
and Visik and Fursikov j.'i2j . 

1.2. Setting of the problem. In the sequel we will restrict to the two dimen- 
sional case, since the nonlinear character of the equations imposes some restriction, 
see Remark 12.51 Furthermore, we fix the values of L and of the Reynolds number to 2 
due to the fact that we will not deal with the singular limit Re — > oo (regarding this 
limit see also the recent works of Clopeau et al. [7J and Mucha [221 ) - In our case we 
will study the following boundary-initial value problem: 

d t u- V -D(u) + (u- X7)u + X7p = f infix (0,T) 

V • u = in fi x (0, T) 

u-n = 5 a g{x,t) onTx(0,T) (1.6) 

u-T + Sn-D(u) -t = onrx(0,T) 

u(x,0) = uq(x) in fi. 
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Note that a similar problem, involving the Smagorinsky-Ladyzhenskaya turbulence 
model together with a nonlinear dependence of (3 on u, has been studied for instance 
by Pares |2S], but in this reference the normal datum g is not allowed to depend on 
the time variable. In the sequel our main interest will be to find weak hypotheses 
on g(x, t) with respect to the time variable (without any essential restriction on the 
space regularity) that allow to prove existence of weak solutions to the Navier-Stokes 
equations, see Theorem ll.il In particular, in our analysis we will focus on two main 
points: 

a) to show the existence of weak solutions in the sense of Leray-Hopf (since we 
do not want to deal with any weaker concept of solution) ; 

b) to use only elementary tools of functional analysis, in order to keep the paper 
intelligible to non specialists. 

In other words, we want to consider solutions in an usual sense and we want also to 
interact with applied people interested in this problem, still keeping all the mathe- 
matical rigor needed to deal with the problem and a certain sharpness of the results. 
In the case of non-homogeneous no-slip conditions, several results of existence and 
uniqueness of suitable solutions can be found in Amann 0. 

Our intent to have a non- vanishing normal datum can be heuristically understood 
also with the following argument: suppose that (for simplicity in two dimensions) we 
have a fictitious boundary Ti and we want to impose a condition on it in order to 
resolve numerically the equation in a smaller domain f2i C f2, that rules out the 
boundary layer (see the figure below). 




Boundary Layer 



Boundary Layer 



r 



B 



Fig. 1. The fictitious boundary 

We have to require, by the incompressibility of the flow, that 
/ W-udx= u-nda = Q, 

J{ABCD} J d{ABCD} 

for each (also curvilinear or infinitesimal) "rectangle" {ABCD} touching the bound- 
ary r as in the figure. Since the behaviour of the flow is not known, in general we 
have 3 



u ■ n ds 



{CD} 



u ■ nds - 



{BC} 



u ■ nds 



{DA} 



?0. 



This may justify the introduction of a non vanishing normal flux, also with very low 
regularity properties, namely the same shared by the trace of a turbulent flow in the 
boundary layer region. 



3 Note that the line integral over the segment {AB} vanishes, since on the "true boundary" V 
both the Navier or no-slip conditions impose that u ■ n = 0. 
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1.3. Main results. In this section we briefly enunciate the results we will prove, 
together with their precise and rigorous statement. 

In the sequel will a bounded, connected, open set in R 3 , locally situated on one 
side of its boundary T, a manifold of (at least) class C 1,1 (Lipschitz-continuous first 
derivatives). The existence of the unit outward normal n derives by results proved in 
Necas 

The first result we will prove is an existence and uniqueness theorem for weak 
solutions of the Navier-Stokes system l)1.6|) . with boundary conditions 1)1. 4|) . 

Theorem 1.1. Let be given g € H 2+ e (0, T; H 3 (r)), satisfying the compatibility 
condition (|1.5|) . / <G i 2 ((0,T) x fi), and assume that uq 6 L 2 (fl), with V • uq = 0. 
Then, there exists a unique weak solution 



Next, we want to study the behaviour of the solution to problem (|l.ti|) as the 
small parameter S converges to zero. In view of the considerations of the previous 
section, one can expect that, as the boundary layer becomes thinner and thinner, the 
solutions will look more and more like the classical solutions corresponding to the 
no-slip boundary condition. Indeed, this is the case, as it is shown by the following 
theorem. 

Let us be the solution of 1)1.6)) (we emphasize the dependence on 5 in this frame- 
work) and let v be the solution to the Navier-Stokes equations with the same initial 
value and no-slip boundary conditions (to be more precise, the vector field v is the 
solution to system 1)3. 2J) )■ As we shall see in Section 



so that the "no-slip solution" represents the average behaviour, once one neglects the 
effect at the boundary. The term us — v can be seen as the "fluctuation term" , which 
takes into account the non-trivial dynamics at the boundary. 



Theorem 1.2. Assume u £ H 1 ^), with V • u = 0, g S H^ +e (Q,T; (T)) 
satisfying the compatibility condition p. 5(1 . and f G L 2 ((0,T) x fl). Then 



In particular, us converges to v in I/°°(0, T; L 2 (fl)) and L 2 (0,T; _ff 1 (f2)). 

2. A result of existence and uniqueness of weak solutions. In this section 
we prove Theorem ll.il For the sake of simplicity, we consider the normal datum as 



namely, dropping the dependence on 6, since it is not relevant in view of the existence 
and uniqueness result we are going to show. In the last section we shall see how to 
deal with a right-hand side that scales by a power of S. 
Let us consider the evolution problem 



u S L°°(0, T; L 2 (fl)) n L 2 (0, T; i? 1 (f2)), 



of problem 1)1. 6)1 . 



u s = v + (D(S3), 




u ■ n = g(x, t) 



' d t u- V • D(u) + (u 
V ■ u = 
< u ■ n = g{x, t) 
5 n • D(it) • t + u ■ t 
u(x, 0) = uq(x) 



V) u + X7p = f in O x (0, T) 



on T x (0, T) 

in f2, 



in SI x (0, T) 
on T x (0, T) 



(2-1) 
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with / £ L 2 ((0,T) x SI) and V • / = 0, just to avoid technicalities, and with g not 
very smooth, say 

gG H^ 2+e (0,T;H^ 2 (T)) (2.2) 

satisfying the compatibility condition 11.511 . 

2.1. Function spaces. We use the classical Lebesgue spaces and in particular, 
we will work exclusively in the Hilbert framework, so we will use the space L 2 . For 
simplicity we do not distinguish between space of scalar, vector, or either tensor valued 
functions, and the symbol | . | will denote the norm in L 2 (Q). The norm in L 2 (T) will 
be denoted by || . ||r- 

In the sequel we will use the customary Sobolev spaces for which we refer to 
Adams £Q and we will use the notion of "trace" over the boundary T of f2. Mainly we 
will use the space iJ 1 (fi) with norm denoted by || . ||#i and its trace space 7? 1 / 2 (r), 
with norm || . ||i/2,r- 

In addition, we define the spaces 

H = {u e L 2 (n) I V -u = 0, u-n = 0onr}, 

and 

V = {ueH 1 {Q)\V -u = 0, u-n = 0onr}, 

and we endow V with the norm \\u\\v = ||Vu||. Moreover, we define the space of 
tangential vector fields as 

ffj={ue H 1 (Q)\u-n = on T}. 

2.1.1. Fractional derivative. In order to define properly the spaces we will 
use, we need also to define fractional derivatives. The fractional derivative may be 
defined through singular integrals 



* U(s,x) , , „ 

-, r — as for < a < 1 , 

(t -a)" ~ 



but for our purposes it is better to deal with a definition via the Fourier transform. 

Given 4>(x,t), defined for t £ [0, T], with values in the Hilbert space X and 
integrable (in the Bochner sense), we define 

*{t, X ) = {^ X) i0rte[ °' T] (2-3) 

I elsewhere, 

and its Fourier transform (with respect to the time variable) is 

#c,0 = --L / 4>{x, t)e- u Ut, 
V Ztt J k 

so that we can define the fractional Sobolev spaces of functions having a— order deriva- 
tive in L 2 : 

iT(lrl;X) := {/ £ L 2 (R;X) : J |C| 2q ||7(0II| « < «)}. 
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2.2. The linear stationary problem. The first step to solve ()2. ID is to consider 
the linear stationary problem 

'-V-D(G)+Vn = inf2x(0,T) 

V • G = in SI X (0, T) . iS 

x ' (2.4) 
G-n = g(x,t) onrx(0,T) v ; 

Jn-~B(G)-T + G-T = onrx(0,T), 

where the time variable is now just a parameter. 

Theorem 2.1. Let be given g G H 1 / 2+e (Q,T;H 1 / 2 (Y)), satisfying the compatibil- 
ity condition (|1.5fl . Then, there exists a unique G solution of (|2.4H such that 

G(x, t) e 1/2+e (0, T; H 1 ^)). (2.5) 

Moreover, there is a constant Co, depending only on SI, such that 

l|VG|| + ||n|| <Co(l + 5-i)\\g\\ Hi{ry (2.6) 

Proof. See Solonnikov and Scadilov and Beirao da Veiga pj- In fact, for 
each i it is possible to solve a linear stationary Stokes problem (with the appropriate 
boundary conditions) that has a unique solution belonging to H (SI). The regularity 
in the time variable is inherited by the function G. 

We give a formal (but completely justified) argument for the estimate (|2.fill . fol- 
lowing the approach to the existence in L 2 -spaces introduced in Beirao da Veiga (SJ 
to find appropriate estimates on G. We consider the bilinear form 

B(u,4>) = I B(u)B(4>) dx 
Jn 

and the functions (6?, II) , that solve l|2.4(l , must satisfy 

B(G,(f>)- [ UV-(j)dx+-- [ G-cpda = 0, Vcj) € Hi. 
Jn o Jr 

In order to deal with the inhomogeneous problem, we introduce a vector field Gi such 
that 

'V-Gi=0 in SI, 

G\ ■ n = g on T, 

\\G^<C\\g\\ Hi{r) . 

The construction of such a vector field is rather standard and can be found for instance 
in Galdi [13]. 

By defining G — G\ + G2, the function G2 must satisfy 

B(G 2 ,(t>)- I nV-(/>dx+- f G 2 -4>da = -B(G x ,<t>) -- f G x -<j)da, 

for each </> tangential to the boundary. If <fi = G2, we get (since Gi • n = and 
V-G 2 = 0) 

l|VG 2 || 2 + i||G 2 || 2 < ||VG 2 || HVGiH + i||G 1 || r ||G 2 ||r, 
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and consequently 



VG 2 || 2 



26 



\G 



2 M r 



< 



IVGxI 



26 



This finally implies that 



|VG 2 || 2 <G[ 1 



'i?2(r) 



where the constant C depends on Q, but it is independent of 6. The estimate on the 
pressure can be obtained by approximation, by studying a slightly different equation 
(see jS] for details). □ 

Indeed, notice that, for our aims, it is enough to have 



G(x, t) € H^ 2+e (0, T; L 2 {Q)) n L 2 (0, T; H\Q)). 



(2.7) 



but currently we do not know the minimal assumption on g in order to have the above 
regularity. Regarding the usual no-slip boundary conditions see the result proved by 
Fursikov, Gunzburger, and Hou |12) . 

2.3. The linear evolution problem. The next step for the analysis of the 
nonlinear evolution problem 1|2.1[) is the following linear evolution problem 

'd t z- V-D(z) + Vg = inOx(0,T) 

V-z = infix(0,T) 

z ■ n = g on T x (0, T) (2.8) 

5r ■ T>(z) ■ t + z ■ t = onrx(0,T) 

z{x, 0) = G{x, 0) inO. 

We shall treat the non-linear problem as a perturbation of such a linear system. Let 
us introduce the new unknowns 

Z(x,t) — z(x,t) — G(x,t) and Q(x,t) = q(x, t) — H(x,t) 

so that we are reduced to a homogeneous problem for the new unknowns (Z, Q): 



d t Z - V-D(Z) + Vg = -d t G inf2x(0,T) 

V-Z = infix(0,T) 

Z-n = onrx(0,T) (2.9) 

5n-D(Z)-T + Z-T = onrx(0,T) 

Z(x,0)=Q mil, 



The above problem is not completely standard, since the right-hand side does not 
satisfy the usual properties. For instance, one can note that dtG does not belong to 
the domain of the Stokes operator, since dtG ■ n ^ 0. This is the main difficulty: the 
low regularity of this term can be treated in a more standard way, while the above 
fact is responsible for a different approach. 

Theorem 2.2. Assume that (G, LT) is a solution to system (|2.4|l . with G satisfying 
the regularity property l|2.7ll . Then, there exists a unique solution (z,q) (where q is 
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unique up to an additive constant not depending on the space variable) to system (|2.8[1 
such that 

z 6 L oo (0,T;I 2 (O))nL 2 (0,T;ff 1 (fl)) nH?- s (0,T;L 2 (n)). (2.10) 



Moreover, 



sup \\z(t)\\ 

0<t<T 



mZ )(t)\\ 2 + l\\z(t)-T\\l)dt+\\z\\^ 



5~ e (0,T;L 2 (n)) 
T \ 



< 



ffS +e (0,T;L 2 (J2)) 



|D(G)|| 



(2.11) 



Proof. By virtue of Theorem 12.11 it is enough to prove the same claim of this 
theorem on the solution (Z,Q) of problem (|2.9|) . Since we just know that dtG 6 
H-i+ E (0,T; L 2 (J1)), we introduce a sequence G N G -ff 1 (R; L 2 (il)) of approximate 
functions such that 

(a) G N \ [0T] — > G in H? +s (0, T; L 2 (£l)), as -* oo, 

f&j ||ftG^||i»(o,r;i»(n)) =JV- 
The way to do this extension is rather standard: first we can define G : IR — > L 2 (f2) 
with an extension by reflection. Then, we consider a sequence /9/v of mollifiers and 
the function Gn will be the restriction on [0, T] of the function p^ * G. 

The proof is based on the Faedo-Galerkin procedure. By Clopeau et al. 0, (but 
also the recent abstract results in [3]) we know that there exists a basis (0 n )neN of 
functions in iJ 3 (SI) of the space V (and also of H), such that 

5n-D(<j) n ) - T + cj) n -T = 0. 

Now, let Zn(t,x) = Y^k=i (n kitf&kix) be the solution of the following (finite dimen- 
sional) linear system of Ordinary Differential Equations (ODE): 

f M) + (BiZnWM) + ]J r (Z^r) (fa-T)da = ~j t {G N ,<f> k ) 

for t € (0,T) and k = l,...,n, where (•,•) denotes the usual L 2 -scalar product. 
Notice that the divergence-free constraint and the boundary conditions on Z^ are 
automatically verified. By using a standard argument it is easy to prove that such 
a system of ODE has a unique solution Z% <E L°°(0, T; L 2 (fl)) n L 2 (0, T; H 1 ^)): 
indeed, by multiplying each equation by the corresponding term Cnk> summing over 
k and integrating by parts over Q, one easily obtains the following estimate 

sup ||^(t)|| 2 + / T (||D(^)|| a + i f T \\Z^\\l)dt<C\\G N \\^ T . L , m , 

0<t<T Jo V "JO J 

for a constant C, depending only on f2. Unfortunately, such an estimate, beside 
being uniform in n, is not uniform in N due to property (b) of the approximate 
sequence (G n )ng~n- Hence, we need other a-priori estimates on the solutions Z^ of 
the finite dimensional problem. We continue working on the Zj? , since such functions 
are smooth enough for the computation that will be performed. 
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We multiply again the equations by the terms Cnk^ sum over k and integrate by 
parts, but we estimate the right-hand side in the following way: 

1 



sup \\Z»(t)\\ 2 - 

0<t<T Jo 



N\\\2 



< \\d t G N \ 



n 7JV11 



Jo 



d t G n 



7N 



dxdt 



H-3+ e (0,T;L 2 (n)) n ™ u H?- e (0,T;L 2 (n)) 



< \\G 



N 



\z 



N 



(2.12) 



so that we only need to show a uniform estimate (with respect to both n and N) of Z^ 
in the space H2- e {Q,T-L 2 (n)). We shall use the Fourier transform characterization 
of the norm of fractional Sobolev spaces (see Adams pQ) to get such estimate. Each 
Z^ (such functions have been defined in (|2.3|0 is a solution of the following equation 



1'^ 



rN 



T>(Z%) • D( 



; z 



N 



dt 



G ■ 4>k + 5{t) G N (0)-fa-S(t-T) {Z»{T) + G»{T))-<j> k , 



for each k = 1, . . . , n, in the sense of distributions with respect to the time variable. 
Here 5(-) is the usual Dirac's delta function. In the frequency Fourier variable £, the 
above equation reads: 



Jn Jn Jr 

= i£ / G N ■ fa + [ G» • fa - e-* T [ (Z»(T) + G N (T)) ■ fa, 
Jn Jn Jn 

see for instance Lions [20], where this tool is used to prove estimates on the fractional 
derivative of the solution. Note that in that reference, and in all involving fractional 
derivatives for the Navier-Stokes equations, the starting point is the existence of a 
weak solution and on it it is possible to prove additional estimates. In our case the 
existence of weak solution derives from the fractional derivative estimates and it seems 
not possible to prove the usual existence results. 
Consequently, we get 



- i£||^(£)|| 2 + ||D(^)(£)|| 2 + -11^(011? = 

= it[d N -ZK+ [ G N (0)-Z%-e-^ T [ (Z^(T) + G N (T))-Z^. 
Jn Jn Jn 



Take the imaginary part and multiply both sides of the previous formula by |£| 2 
with a < I so that, by using Young's inequality, one gets 

|£Hl^(£)H 2 < C\e a \\G N \\ 2 + C\e a - 2 (\\G N (T)\\ + \\Z?{T)\\ + \\G N (Q)\\f 
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In order to estimate the integral J R |£| 2a ||Z^(£)|| 2 we split it in two parts; by the 



above estimate, 



\cn\z^o\\ 2 < 

ci>i 

<C [ lei^llG^f + CdlG^^II + II^^H + IIG^WII) 2 / |C| 2Q - 2 ; 
the first term on the right-hand side is controlled by Gil G w II 2 i , , while 

6 J " "ffI + t (0,T;I ! (!l)) 

||Z^(T)|| 2 < C\\G \\ H i+. {0T . L2{n)) \\Zn\\ H i-e i0T . L2(n) y 

by virtue of 1(213)1 : 11(2^(0)11 is bounded by ||G(0)||, and finally, by using the Morrey 
inequality that implies iJ 1 / 2+e (0,T) C G([0,T]) (see Adams jf]) we get 

\\G N (T)\\ < \\G N \\ H ^ (0 T . L2(n)y 

The second part is estimated as follows, by using Parseval's theorem, Poincare in- 
equality and estimate (|2.12() . 



/ \e a \\zz\\ 2 dt< [ \\z»fdt= [ wz^tWdt 

J\£\<1 JH -JO 

<C T \\Y>{Z»)fdt 
Jo 

- H ff 3+ E (0,T;L2(f2))ll^n Hfl-3- e (0,T;L2(n))- 

In conclusion, by collecting all of the above estimates, we finally get that, for each 
e G (0, |), there exists a constant G, depending only on f2 and e, such that 

\\ Z n \\ H h-*(0,T;LZ(Cl)) " C|1 G " 11 ff* +' (0,T;i=(O)) ' (2 ' 13) 

that, together with 12.12|l . says that is bounded, uniformly in n and N, in the 
spaces Hi~ e {0,T;L 2 {VL)), L°°{Q, T; L 2 (Q)) and L 2 (0, T; i? 1 (SI)). 

Hence, it is possible to extract a (diagonal) subsequence converging weakly in 
L 2 (0, T; i? 1 ^)), weakly-* in L°°(0, T; £ 2 (S1)) and strongly in L 2 ((0,T) x SI) to the 
unique solution Z of problem [|315)l . Indeed, Z G iJ 2~ £ (0, T; L 2 (S1)), that is, the 
topological dual space of ff~2+ E (0, T; L 2 (S1)), the space to which d t G belongs to. 
Furthermore, by passing to the limit and by using the semi-continuity of the norms, 
it follows that Z satisfies the claim stated at the beginning of this proof. □ 

Remark 2.3. The assumption g G H 2 +e (0, T; H 2 (r)), which in turn gives 
G G H 2+^(0, T;H\n)), seems to be rather technical, for the presence of e. If 
g G Hi (0, +00; Hi (T)), it follows that G G #5(0, +oo;H 1 (n)) and Theorem BJ 
holds accordingly. 

Indeed, the main point is estimate l|2.12|l . in which the right-hand side becomes 
\\G N \\ 1 ||Z„|| 1 and, following the lines of the proof presented above, the esti- 
mate on the Fourier transform gives that Z^ is bounded, uniformly in n and N, 
in H 2 (0, +00; L 2 (fl)). Notice that in the critical case H^, we work on the whole time 
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interval [0,+oo), to avoid the boundary terms \\G N (0)\\ and \\G N (T)\\, that cannot be 
estimated by using the Morrey inequality. 

We also note that this small relaxation on the assumptions on G requires to add 
the hypothesis G S L°°(0, T; L 2 (fl)). Otherwise the function z will not belong itself 
to L°°(0, T; L 2 (fl)) and this fact is crucial to prove the corresponding bound for weak 
solutions to the full nonlinear Navier-Stokes problem. 

2.4. The non-linear problem. In this section we finally prove Theorem 11.11 
Again, we make use of an auxiliary problem; namely, we introduce the new variables 

U = u — z and P = p — q, 

where (z, q) is the solution to the linear evolution problem (|2.8(l . and the pair (17, P) 
solves the following problem 

d t U- V ■D(U) + [(U + z) ■ V](U + z) + VP = / infix (0,T) 

V • U = in fl x (0, T) 

J7-n = onrx(0,T) (2.14) 

5n-T>{U)-T + U -r = onrx(0,T) 

^ U(x, 0) = u (x) ~ G(x, 0) in fi. 

By virtue of Theorem 12. 21 the existence Theorem 1 1.1 1 for the non- linear problem is a 
straightforward consequence of the following proposition. 

PROPOSITION 2.4. Assume that (G, II) is solution to system (|2.4|) . with G £ 
H?+ e {0,T;L 2 {Q))C]L 2 {Q,T- 1 H 1 {VL)), then there is a unique 

U E L°°(0, T; L 2 {Vl)) n L 2 (0, T; F x (0)) (2.15) 

solution to problem l|2.14|l . Moreover, the following estimate holds true: 

sup \\U(s)\\ 2 + f (}\T>(U)f + h\Uf T )ds< 

0<s<t Jo V J 

< \\u - G(-,0)\\ 2 e A ^ + C [ (||/f + ||Vz( S )|| 2 ||z( S )||)e^- A ^ ds, (2.16) 

Jo 

where 

A(t) =Ct + c(l + ||*||£o. (0 ,T;L»(n))) J W^zfds, 

and C is a constant depending only on Q. 

Proof. The proof is rather standard and proceeds via a Faedo-Galerkin approx- 
imation, as in the proof of Theorem 12.21 We only show an a-priori estimate, whose 
computations are formal, but are completely meaningful at the level of the Faedo- 
Galerkin approximate functions. Multiply l|2.14JI by U and integrate by parts, to 
get 

\j t \\U\\ 2 + ||ID(L0|| 2 + \\\U\\ 2 = jf U-{(U + z) ■ V](U + z) + jT / • U. 

The estimate of the integral involving / is straightforward, since it is bounded by 
Il/H 2 + ||J7|| 2 . We estimate the non-linear term in the right-hand side by using the 
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Gagliardo-Nirenberg inequality 4 , 



IMU« < c\\u\\ 1/2 \\\!u\\ 112 v« e H\n). 



(2.17) 



We first observe that, since V • U = and U ■ n = 0, 




U ■ (U ■ V)U = 



and / 

Jo. 



U-(U-V)z 




z-(U ■ V)Z7, 



so that, using repeatedly the Gagliardo-Nirenberg inequality above given, Holder's 
inequality and Young's inequality, we get 



Since, by Theorem z € L°°(0,T; L 2 (Cl)) n L 2 (0, T; both terms (1 + 

||z|| 2 )| Vz|| 2 and ||Vz|| 2 ||z| are integrable in time and, by Gronwall's lemma, we 
can deduce that U is bounded in L°°(0, T; L 2 {VL)) and i 2 (0, T; H 1 ^)). Moreover, 
formula l|2.16|l also follows. 

Finally, uniqueness of the solution follows from similar arguments. Indeed, if U 
is the difference between two solutions U\ and f/a, one easily gets 



and, since ?7(0) = 0, from Gronwall's lemma it follows that U = 0. □ 

Remark 2.5. In the proof of the result of this section we used in a fundamental 
way estimate (I2.17|) . In the three-dimensional case this inequality is not anymore 
true. Instead, it holds 



and the latter can be used to prove just local existence of weak solutions. The global 
result proved in the 2D case stems in an essential manner on the stronger estimate 
and this is the critical difference between the two cases. 

3. The proof of Theorem ll.21 In this section we study the limit of the solution 
of the vorticity seeding model as 5 — > 0. We prove the convergence of the solutions 
of to the solutions of the non-stationary Navier-Stokes system with the no-slip 
boundary condition (|3.2() . This supports the idea of using a non-standard boundary 
condition on the new boundary T% C ft, such that the region between T and Tx is very 
narrow. When the width of this region (presumably the boundary layer) shrinks to 
zero, then the solution of the Navier-Stokes equations with the usual no-slip boundary 
condition is recovered. 

4 Note that such a inequality is a little bit more general than the so-called Ladyzhenskaya in- 
equality, since the functions are not vanishing on the boundary of f!. 




<-||vc/|| 2 + q|z||l 4 ||c/|| 2 + q|z||| 4 ||Vz||t||c/||l 

< ^l|VC/|| 2 + C||z|| 2 ||Vz|| 2 ||?7|| 2 + C||z||i||Vz|| 2 ||[/||i 

< \\m\\ 2 + C(l + ||z|| 2 )||Vz|| 2 n| 2 + C||Vz|| 2 ||z||. 



^||C/|| 2 <(||C/ 2 ||| 4 + ||z||l 4 + ||Vz|| 2 )||C/|| 2 
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3.1. Comparison of solutions with different boundary data. Let g S 

Hi+ e (0,T;Hi(T)), satisfying the compatibility condition ifO))) . / E L 2 ((0,T) x fi) 
and uo £ H 1 ^), with V ■ ito = 0. Without loss of generality we can assume that 
uq = 0. Denote by Gs the solution of the linear stationary problem (|2.4fl . with 
boundary condition • n = S g. Consider the solution (us,ps) to the system 



d t u s - V ■ D(u a ) + • Vu 5 + V Ps = f in fi x (0, T) 

V • u s = in fi x (0, T) 

u 5 -n = 5g(x,t) onT x (0,T) 

5n-7D(us)-T + u s -T = onTx(0,T) 

us(x, 0) = uo(a;) + G(x, 0) in fi, 



(3.1) 



We want to show the convergence of the vector valued function us to the solution v 
of the Navier-Stokes equations with zero Dirichlet data: 



'd t v- V • D(v) +v ■ Vw + Vtt = / infix(0,T) 

V-i> = infix (0,T) 

v = on T x (0, T) 

t v(x, 0) = uq in fi. 



(3.2) 



To this end, we also introduce the solution zs to the linear evolution problem H2.8(l . 
with boundary condition zs • n = S g, and we set 

Us = us — zs and ws = us — z$ — v. 

The function ws satisfies the following homogeneous system 

~d t w 5 - V ■B(w s )+R(w 5 ,z s ,v,U s )+Vr = infix (0,T) 

V- 10,5 = infix (0,T) 

w 5 -n = onrx(0,T) (3.3) 

5n • B(ws) - t + ws • t = -<5n • D(«) • r on T x (0, T) 

ws(x, 0) = in fi, 

where 

R(w s , z s , v, Us) = (Us ■ V)w s + (w 5 ■ V)« + (U s ■ V)z s + (z s ■ V)U S + (z s ■ V)z 5 . 
We multiply the first equation in l|3.3|) by ws and integrate by parts to get 



Id 

2d7 KI 



\D(ws)\\ 2 + ^\\w s \\l = - J ' (n-D(v)-T)(w s -T)da- J w s - R{w s , Z S ,v,U s ). 



The boundary integral in the right-hand side may be increased in the following way 



J (fi ' • z)(ws ■ L) da 



< ^\Hh + ^IKIIr- 



(3.4) 



We analyse then the second integral. By integration by parts over fi (note that 
Us ■ n = on T) and by using Holder's inequality, the Gagliardo-Nirenberg inequality 
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and Young's inequality, we obtain 



w s ■ R(ws, z$, v, Us) dx 



< 



< i||D( W5 )|| 2 + C[||V«|| 2 + ||VC/ 5 ||I||Vz d -||t + \\WUs\\^\\Wz s \\i + \\Wz s \\ 2 ]\\ws\ 
+C[||Vt7 4 ||'||Vz a ||'||£7 4 || + ||VI/5||*||VZ5||»|M| + ||Vz 5 || 2 ||z 4 ||]. 
Indeed, for example, 



ws ■ (z s ■ V)Us 



< IKIU*||*«IU*||vi7|| 

<C\\ws\\^\\T>(ws)\\^\\zs\\^\\D{zs)\\^\\\7Us\\ 

< ^\\lD(ws)\\ 2 + C\\VUs\\?\\Vz s \\i\\ws\\l\\zs\\i 

< h\D(ws)\\ a + C\\VU B \\ * || V^|| § (IKI| 2 + \\z s \\). 

o 

For simplicity we set 

ij>{t) := C[\\VUsf*\\Vzstz\\Us\\ + \\VU 5 \\^\\Vz s \\^\\z s \\ + \\ Vz s \\ 2 \\zs\\] , 
<f>(t) := [||Vi;|| 2 + ||Vl7a||f ||V*«||* + \\ VU S \\ > || Vz 5 \\ § + ||V2 5 || 2 ], 

x(t) HM&», 

so that, by recollecting all together the estimates we have before obtained, we end up 
with the following differential inequality 

^|K|| 2 + ||JDK)|| 2 + i|K||J < 6 X (t) + m + <P(t)\\ws\\ 2 (3.5) 

3.1.1. Estimate of t/j and <f> in terms of 8. First, we note that the functions 
ip and <fi belong to L 1 (0,T), by virtue of the regularity properties (|2.1UI) and l|2.15[) of 
zs and Us, respectively. The function %(t) belongs to i 1 (0,T) as well, since 

\\v\\mo,T:H2) < C[\\V(X,0)\\ H 1 + ||/||z, 2 ((0,T)xfi)], 

see for instance Constantin and Foias, 9 . And also in this point is crucial to consider 
the 2D problem, since such estimate is available just for small times in the 3D case. 

It is important now to sharply check the correct behaviour of the functions <p and 
in terms of S. Before going further, we collect in the following lemma the results 
of Section [21 that we shall need. 

Lemma 3.1. Under the assumptions of Theorem \l.Hl the following estimates hold 

2 i / IItp>^-M|2 J + n Xll„l|2 



(i) sup ||z,|| 2 + / ||B(z 5 )||^t<C 1( %| ■ , 

0<t<T Jo H ( '^ fl5 (')) 

r T i 

(ii) sup \\Us\\ 2 + / (\mUs)\\ 2 + -JUs\\ 2 r )dt<C 2 , 

0<t<T JO K 6 / 

for all < S < 1, where the constant C\ depends only on Q and the constant C : 
depends only on fi, T, ||/||l2, and \\g\\ i +e . 
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Proof. Property (i) is a consequence of Q2.1ip . Indeed, by (|2.6|1 and Poincare 
inequality, 



< C\\B(G B )\\ < CV5\\g\\, 

Again, from l|2.6|> . 



so that \\G s \\ Hi+ a)) < C->/S||s|| ff4+ . T;J . 2(n)) . 



\\T>(G 5 )\\ 2 dt < C6\\g\\ 2 1 + 

and in conclusion (%) holds true. 

As it concerns (ii), we have that, by Poincare inequality, i|2.(j[l and Sobolev em- 
beddings, 

||G 5 (-,0)|| < C||E>(G?,)(-,0)|| < CV5\\g(;0)\\ < CV~5\\g\\ H ^ (nT . L2 y 
Moreover, since 5 < 1 and from the above estimates it follows 

A(t)<A(T)<CT + C(l + \\g\\ 4 1 , ) 

that is, A(t) is uniformly bounded by a constant independent of <5, and 

\\\ff + ||Vz 5 || 2 ||z 5 ||)e^)-^) ds < e^( f Il/H 2 + ||z 5 || iO o (0iT;L2) ||D(z 5 )|| 2 ) 
o Jo 



<(ll/ll^((o,T )x n) + C|bll^ + Je A(T) 



as well, so that, by H2.16JI . also (ii) follows. □ 

The first consequence of the above lemma is that 



f <t>{s) ds < f || V«|| 2 dt + C(fi, T, /, g), 
Jo Jo 



with a bound uniform in 5, for S small. Moreover 

I i 
^(s)ds<C[ [ \\VUs\\ 2 ) If ||Vz 5 || 2 ) ||^|| L o 



2 1 \ 3 / .y 

2 I / / ||V7~ l|2 



+ C^y ||V^|| 2 J \\Vz s rj \\z s \\l°o +C\\z s \\ l ~ j \\Y:,\ 

<C(Si+Si+Si) 
<C6^, 

so that, by l|3.5|) and Gronwall's inequality, it follows that 



sup 

0<t<T 



|M| 2 + f T (||DK)|| 2 + JlMI 2 ) dt < [ T [5 X (s) + i;(s)}e^^ dr ds 



< e I Q T <!>(*) ds / [6 x ( 8 ) + il>(8))]d8 

Jo 

< C8i 
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and finally 

sup |M| 2 + [ T (\\B(w s )f + h\w s f r )dt = 0(S^. 

0<t<T Jo v J 

Since the same norms of z$ are of order 8 i, by the previous lemma, we finally deduce 
that ug — v is of the order 8$ , and the theorem is proved. 

Remark 3.2. The tangential trace of the function w$ converges a little better, 
since from the above estimate we can deduce 

IN|| L 2 ( rx(o,T)) =0(S 5 / 6 ). 

Remark 3.3. The estimates show that there is a crucial loss in the estimates 
due 1) to the boundary effect, in the estimate Q2.6J1 : 2) to the non-linearity (recall the 
contribution of Us in the estimate of ip(t)). We observe that the original model was 
supplemented by the following boundary condition 

u ■ n = S 2 g, 

but we used in fact u ■ n = 8g. By running through the proof of Theorem 11.21 one 
gets, using the boundary condition 

u ■ n = 8 a g, 

that zs = 0(8 a -i), so that fip(t)dt = 0(<0(*-*)) and finally w s = 0{8^ a ~^). 
In particular, for a = 2, that is, the value corresponding to the original model, the 
behaviour of zg matches the loss in the estimates due to the boundary term Ij3.4|l . 
Hence, for larger values of a, the convergence remains slow, because of this term. 
Supposedly, the loss in the convergence rate caused by this term is an intrinsic feature 
of the problem, while those caused by the non-linearity may have technical reasons. 
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